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Abstract 

We study thermal leptogenesis in classes of neutrino mass models based on the 
seesaw mechanism with three right-handed neutrinos and sequential right-handed 
neutrino dominance. The flavour-dependent Boltzmann equations are solved ap- 
propriate to both the Standard Model and the Minimal Supersymmetric Standard 
Model. Within these classes of models we investigate constraints and expectations 
on the individual decay asymmetries and washout parameters from the present 
data on neutrino masses and mixings. In many cases of physical interest flavour 
effects are shown to have important consequences for the estimation of the pro- 
duced baryon asymmetry in leptogenesis. We also establish and analyze the link 
between the leptonic CP violating phase 6, observable in neutrino oscillations, and 
the CP violation required for leptogenesis, where flavour-dependent effects have a 
significant effect. In general our results show that flavour-dependent effects cannot 
be ignored when dealing with three right-handed neutrino models. 
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1 Introduction 



Thermal leptogenesis [Ij is an attractive and minimal mechanism to generate the baryon 
asymmetry of the Universe n-B/n^ ~ (6.10 ± 0.21) x 10~^°, or, normalised to the entropy 
density, Yb ~ (0.87 ± 0.03) x 10~^° [2j. A lepton asymmetry is dynamically gener- 
ated and then converted into a baryon asymmetry due to {B + L)-violating sphaleron 
interactions [3J which exist in the Standard Model (SM) and its minimal supersymmet- 
ric extension, the MSSM. Leptogenesis can be implemented within the "Seesaw" (type 
I) model [1], consisting of the SM (MSSM) plus three right-handed (RH) Majorana 
neutrinos (and their superpartners) with a hierarchical spectrum. In the simplest case, 
the lightest of the RH neutrinos is produced by thermal scattering after inflation, and 
subsequently decays out-of-equilibrium in a lepton number and CP-violating way, thus 
satisfying Sakharov's constraints [3]. 

The asymmetry is commonly calculated by solving a Boltzmann equation for the 
total lepton asymmetry (flavour-independent approximation) and for the abundance of 
the lightest RH neutrino P, d E] . However, this flavour-independent treatment is rigor- 
ously correct only when the interactions mediated by charged lepton Yukawa couplings 
are out of equilibrium. The impact of flavour in thermal leptogenesis has been first 
addressed in Ref. [S] and then studied in detail [101 [HI [12] , including the quantum oscil- 
lations/correlations of the asymmetries in lepton flavour space [TU]. It was shown that 
the Boltzmann equations describing the asymmetries in flavour space have additional 
terms which can significantly affect the result for the final baryon asymmetry. This is 
because leptogenesis involves the production and destruction of right-handed neutrinos, 
and of a lepton asymmetry that is distributed among distinguishable flavours. The pro- 
cesses which wash out lepton number are flavour- dependent, e.g., the inverse decays from 
electrons can destroy the lepton asymmetry carried by, and only by, the electrons. The 
asymmetries in each flavour are therefore washed out differently, and will appear with 
different weights in the final formula for the baryon asymmetry. This is physically in- 
equivalent to the treatment of washout in the flavour-independent approximation, where 
indistinguishable leptons propagate between decays and inverse decays, so inverse decays 
from all flavours are taken to wash out asymmetries in any flavour. Flavour-dependent 
effects in leptogenesis have been studied in the two right-handed neutrino model |12j . 
where they have been shown to be relevant. It remains to be seen how important these 
effects are in more general neutrino mass models with three right-handed neutrinos. 

The latest experimental data on neutrino oscillations is consistent with (approxi- 
mate) tri-bimaximal mixing [13]. The fact that the tri-bimaximal neutrino mixing ma- 
trix involves square roots of simple ratios motivates models in which the mixing angles 
are independent of the mass eigenvalues. One such class of models are seesaw models 
with sequential dominance (SD) of right-handed neutrinos [Hj. In SD, a neutrino mass 
hierarchy is shown to result from having one of the right-handed neutrinos give the 
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dominant contribution to the seesaw mechanism, while a second right-handed neutrino 
gives the leading sub-dominant contribution, leading to a neutrino mass matrix with 
naturally small determinant. In a basis where the right-handed neutrino mass matrix 
is diagonal, the atmospheric and solar neutrino mixing angles are determined in terms 
of ratios of Yukawa couplings involving the dominant and sub dominant right-handed 
neutrinos, respectively. If these Yukawa couplings are simply related in some way, then 
it is possible for simple neutrino mixing angle relations, such as appear in tri-bimaximal 
neutrino mixing, to emerge in a simple and natural way, independently of the neutrino 
mass eigenvalues. If SD is combined with vacuum alignment in flavour models then 
tri-bimaximal neutrino mixing can be readily achieved [T7]. The SD mechanism has 
been widely applied to a large variety of unified flavour models involving three right- 
handed neutrinos [T3] , and it is therefore of interest to see how important the effects of 
flavour- dependent leptogenesis are in this more general framework. 

In this paper we study thermal leptogenesis in classes of neutrino mass models based 
on the seesaw mechanism with sequential dominance, taking into account lepton flavour 
in the Boltzmann equations. We generalize the issue of including flavour effects in the 
Boltzmann equations to the supersymmetric case. Within this class of model we inves- 
tigate constraints and expectations on the individual decay asymmetries and washout 
parameters from the present data on neutrino masses and mixings. Flavour effects are 
shown to have important consequences for the estimation of the produced baryon asym- 
metry in this class of models. Flavour-independent leptogenesis has previously been 
considered for sequential dominance models in [TSj. The present analysis clearly goes 
well beyond the previous analysis by considering the important flavour-dependent ef- 
fects. Also the link between the leptogenesis phase and the MNS neutrino oscillation 
phase has been explored in sequential dominance in [19]. Here we shall revisit this link 
in the light of flavour effects and obtain new links between the cosmology and neutrino 
oscillation physics. 

The remainder of the paper is set out as follows: In Sec. [21 we review the treatment 
of flavour in the Boltzmann equations and generalize it to the case of the MSSM. Sec. |3] 
contains a review of sequential dominance (SD) in the seesaw mechanism. In Sec. IH 
we investigate leptogenesis in neutrino mass models with SD, taking into account lep- 
ton flavour in the Boltzmann equations. Generic properties of the decay asymmetries 
and washout parameters are derived in Sec. 14.11 and the link between the leptonic CP 
violating phase 6, observable in neutrino oscillations, and the CP violation required 
for leptogenesis is established in Sec. 14.21 Sec. 14.31 contains examples to illustrate our 
results. In Sec. El we conclude. 
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2 Flavour Matters in Leptogenesis 



2.1 Temperatures where Flavour Matters 

In the SM extended by right-handed (singlet) neutrinos Ni {i = 1,2,3) with Majorana 
masses Mj, the additional terms of the Lagrangian are given by 

C = -{X.)a^{r■H)N'-^W{MnRhN'''+ii.c., (1) 

where the dot indicates the SU(2)L-invariant product, ■ H) := £"(ir2)'^^(-ff)6, with ta 
(A G {1, 2, 3}) being the Pauli matrices, ia (a = e, /i, r) are the lepton SU(2)L-doublets 
and H is the Higgs field which develops a vacuum expectation value (vev) of {H'^) = = 
175 GeV in its neutral component after electroweak symmetry breaking. We will work 
in a basis where the charged lepton Yukawa matrix and the mass matrix of the right- 
handed neutrinos are diagonal, i.e. Ae = diag(?/e, y^, 2/r) and Mrr = diag(Mi, M2, M3), 
respectively. We will assume a hierarchical spectrum of right-handed neutrino masses. 
Ml < M2 < M3, in the following. 

In the SM, the fiavour-independent formulae to describe leptogenesis are only ap- 
propriate when the dynamics takes place at temperatures larger than about 10^^ GeV, 
before the charged lepton Yukawa couplings come into equilibrium, estimating the inter- 
action rate for a Yukawa coupling ?/q, as F,- ^ 5 x 10~^ T [20]. However, if leptogenesis 
occurs at smaller temperatures T ~ Mi, where Mi is the mass of the lightest RH neu- 
trino, then one has to distinguish two possible cases. If 10^ GeV ^ Mi ^ 10^ GeV, 
then charged fi and r Yukawa couplings are in thermal equilibrium and all flavours in 
the Boltzmann equations are to be treated separately. For 10^ GeV -C Mi <C 10^^ GeV, 
only the r Yukawa coupling is in equilibrium and is treated separately in the Boltzmann 
equations, while the e and /i flavours are indistinguishable. 

In the MSSM, extended by singlet superfields Np {i = 1,2,3) containing the right- 
handed neutrinos A^* as fermionic components, we use a notation analogous to the SM. 
The additional terms of the superpotential are given by 

W = (A,)„,(^■i^„)iVC^■ + ^iV«(MRR),,iVC^ (2) 

where hats denote superfields. In the MSSM, the vev of the Higgs field H^^, which 
couples to the right-handed neutrinos, is given by = = sin(/3) x 175 GeV, with 
tan P defined as usual as the ratio of the vevs of the Higgs fields which couple to up-type 
quarks (and right-handed neutrinos) and down-type quarks (and charged leptons). 

In the MSSM, the fiavour-independent formulae can only be applied for temperatures 
larger than (1 -|- tan^ (3) x 10^^ GeV, since the squared charged lepton Yukawa couplings 
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in the MSSM are multiplied by this factor. Consequently, charged jj and r lepton 
Yukawa couplings are in thermal equilibrium for (1 + tan^/3) x 10^ GeV <^ Mi <ti 
(1 + tan^/5) X 10^ GeV and all flavours in the Boltzmann equations are to be treated 
separately. For (1 + tan^/?) x 10^ GeV < Mi < (1 + tan^) x 10^^ Qgy, only the r 
Yukawa coupling is in equilibrium and only the r flavour is treated separately in the 
Boltzmann equations, while the e and /i flavours are indistinguishable. 

In what follows we will concern ourselves mainly with the regime where all flavours 
in the Boltzmann equations are to be treated separately, i.e. 10^ GeV <^ Mi <C 10^ GeV 
in the SM and (1 + tan^ f3) x 10^ GeV < Mi < (1 + tan^ /3) x 10^ GeV in the MSSM. We 
will comment on the other two regimes and point out the differences between flavour- 
independent approximation and the flavour- dependent treatment, where lepton flavour 
is taken into account correctly in the Boltzmann equations. We start with the SM and 
then turn to the MSSM. 



2.2 The Boltzmann Equations in the SM 

In the regime where all lepton flavours are to be treated seperately, the Boltzmann 
equations in the SM are given by 



dYj 



— Z 



d^; 



sH{Mi 



-(7d + 7s,Ai=i) 



Yr 



Ni 



VC(i 



- 1 



dz sH{Mi) 



SlA'lD + 7S,AL=l) 



Y, 



\ 2 ~'~ mAL=l) yeq 



(3) 

.(4) 



where there is no sum over a in the last term on the right-side of Eq. (jll) and where 
z = Mi/T with T being the temperature. Yat^ is the density of the lightest right- 
handed neutrino Ni with mass Mi. defined as Y/^^ = Yb/3 — Y^^, where 
Yl^ are the total lepton number densities for the flavours a = e, fi,T and where is 
the total baryon density. It is appropriate to solve the Boltzmann equations for Y/^^ 
instead of for the number densities of the lepton doublets ia, since Aq, = B /3 — La 
is conserved by sphalerons and by the other SM interactions. For all number densities 
Y, normalization to the entropy density s is understood. Y^^ and Y^'^ stand for the 
corresponding equilibrium number densities, is the thermally averaged total decay 
rate of A'^i and 7s,al=i represents the rates for the AL = 1 scattering processes in the 
thermal bath. Notice, in particular, that 7s,al=i contributes to the asymmetry, as was 
recently pointed out in [12]. 

The corresponding flavour-dependent rates for washout processes involving the lep- 
ton flavour a are 7^1 (from inverse decays involving leptons ia) and 7^ al=i- ^'^^ brevity. 
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we have not displayed further contributions from AL = 2 scatterings, which can be ne- 
glected under conditions we will specify below, ei^a is the decay asymmetry of Ni and 
H(T) is the Hubble parameter. 

The matrix A, which appears in the washout term, is defined by = J2f3 ^a/s ^a^- 
The values of its elements depend on which interactions, in addition to the weak and 
strong sphalerons, are in thermal equilibrium at the temperatures where leptogenesis 
takes place. Below 10^ GeV in the SM, A is given by [1^ 



/-151/179 20/179 20/179 \ 
A^^ = 25/358 -344/537 14/537 . (5) 
\ 25/358 14/537 -344/537/ 

Between 10^ and 10^^ GeV in the SM, regarding the leptons only the interaction me- 
diated by the r Yukawa coupling is in equilibrium, and the lepton asymmetries and 
B/3 — La asymmetries in the e and fi flavour can be combined to Y2 = Ye+fj, and 
Ya2 = ^Ae+A^- In this temperature range, A is given by [I2] 

SM _ /-920/589 120/589 



^ \ 30/589 -390/589 ) ' ^^"^ 

Above 10^^ GeV in the SM we recover the flavour-independent treatment, where all 
asymmetries can be combined to Fa = ^Ae+A^+A^; and A is given by A^^ = —1. 

Eqs. and (jlj) can be safely used in the range of temperatures in which the lepton 
Yukawa reactions for each flavour are fully in equilibrium. Indeed, for values of Mi close 
to 10^ GeV in the SM the reactions induced by the muon Yukawa coupling are about 
to be in equilibrium and the quantum oscillations of the asymmetries Fe/x might not 
have been dumped fast enough to be neglected. Eqs. (jlj) may be generalized to include 
quantum oscillations, following Ref. [TU]. However, preliminary numerical simulations 
have shown that the off diagonal terms change the final diagonal lepton asymmetries 
by factors of order unity and therefore, from now on, we will safely neglect them and 
restrict ourselves to Eqs. ([3]) and (HI). 

One can implement the relevant parameters for connecting leptogenesis to neutrino 
models directly in the Boltzmann equations, following [12]. Let us therefore introduce, 
in addition to the decay asymmetries ei^a, the parameters K^, which control the washout 
processes for the asymmetry in an individual lepton flavour a and K, which controls 
the source of RH neutrinos in the thermal bath, as 
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H{Mi) is the Hubble parameter at T = Mi, given by H{Mi) ^ l.QQ^Mf/Mp with 
^SM _ X06.75 being the effective number of degrees of freedom in the SM, X^, denotes 
the neutrino Yukawa matrix (using left-right notation) and T^j^i^ i^Nja^ decay 
rate of Ni into Higgs and lepton doublet £a (or conjugate final states, respectively). The 
thermally averaged decay rate 7/5 is then given in terms of Ftvi^^ by 

id{z) = E , 7S(^) = ^ ^ (r^... + ^N-iJ , (8) 

a 2 V / 

where Ki and K2 are the modified Bessel functions of the second kind. This allows to 
replace 



in Eqs. (E]) and (jl]). Defining in addition two functions /i and by 

+ 7s,AL=i = iDfi , Y + 7w,Ai=i = 7S/2 , (10) 

we can re-write the Boltzmann equations with correct flavour treatment in a simplified 
form as follows [12]: 



The function /i(-2) accounts for the presence of AL = 1 scatterings and f2{z) accounts for 
scatterings in the washout term of the asymmetry [6[ [7] . In our numerical computations 
we only include processes mediated by neutrino and top Yukawa couplings, following 
Ref. [7j. This means, we neglect AL = 1 scatterings involving gauge bosons [221 E] and 
thermal corrections [6], but we take into account corrections from renormalization group 
running between electroweak scale and Mi [9l [23] . We also neglect AL = 2 scatterings, 
which is a good approximation as long as Ka S> 10 x Mi/(10^^ GeV) [12j. Finally, 
according to the usual assumptions for computing the damping rates in the Boltzmann 
equations, i.e. that elastic scattering rates are fast and that the phase space densities for 
both, fermions and scalars, can be approximated as f{Ei,T) = {ni/'rf^)e~^^/'^ ^ where 
= ^TmfK2{mi/T) with gi being the number of degrees of freedom, we use 

vr«^. YS'^A^^/M.). (13) 



6 



The final lepton asymmetry in each flavour is governed by three sets of parameters, 
which can be computed within a neutrino model: ei and K = Ka- e\^a are the 
decay asymmetries of the lightest right-handed neutrino A^^i into Higgs and lepton 
doublet £o, defined as 

^ _ TiVi^. - "^N-i^ ^ ^^^^ 

with the decay rates T^,^^ = r(A^i ^ EJ^) and T^^^^ = r(A^i H*I^). SU(2)l- 
indices in the final state, not displayed explicitly, are summed over. In the SM, the 
tree-level decay rates are 

^ N^U + ^ N{i^ - ^^^1 ^ • y^^) 

The decay asymmetry (arising at one-loop order) is ^211 [25] 

SM_ 1 Ej=2,3l"^[('^t)la[AtA.]lj(Anj»] sMfMj\ . . 

with the loop function g in the SM given by 

'1 + x 



+ 1 - 1 + x In , 

1 — X \ X 



2Vi ■ 



(17) 



Alternatively to and K defined in Eq. ((Tj), the parameters rrii^a and mi will be used 
in following, which we define as 



mi^a = {>^l)la{K)alJ^, mi = ^ mi,o , (18) 

with f u = 175 GeV. They are related to and K by 

= — , or equivalently Ka = — —, with ml^ ^ 1.08 x 10"^ eV . (19) 
m* m* 

2.3 The Boltzmann Equations in the MSSM 

In the MSSM, the density Y^^ of right-handed sneutrinos as well as the densities 
of the slepton doublets have to be included in the Boltzmann equations. Denoting the 
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total (particle and sparticle) i?/3 — asymmetries as ^a^, the simplified Boltzmann 
equations are given by 

^ = -2A-.|||/.(.)(V..-V5^). (20) 



dz 



The matrix A is defined via the relation = X]/3 ^a^.; with = + ^5 being 
the combined densities for lepton and slepton doublets. 

To obtain Eqs. fl2Ul) - fl22|) we have made use of the fact that the tree level decay 
rates satisfy 

with T^-^i^ + Tj^j^ given in Eq. (|T^ . leading to the identities 

i7(Mi) H{M,) H{Mi) H{Mi) H{Mi) H{M{)' ^ ' 

with K, Ka (and mi, mi^^) being defined analogously to the SM case (c.f. Eqs. (J?!) and 
(HE])). From Eq. ([19]), we'find 

m^ssM ^ sin2(/5) X 1.58 x 10-=^ eV , (25) 

using (7^^^^ = 228.75 for computing H{Mi). wi^ggj^ relates K, Ka to fhi, ffiia in 
the MSSM by the analogous of Eq. ( |T9i) . In this conventions, the functions /i and 
/2 in the Boltzmann equations are approximately unchanged (apart from an obvious 
modification in one of the scattering terms, which has only small effects). Note that 
Vu = sin(/?) X 175 GeV in the MSSM. In Eqs. (^01) and (PTI) . processes with particles and 
superpartners are combined, leading to the additional factor of 2 on the right-side of the 
equations. Ignoring supersymmetry breaking, the right-handed neutrinos and sneutrinos 
have equal mass Mi. With the usual approximation of taking Boltzmann statistics for 
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both, fermions and scalars, we have used ^ Y^^ to combine the washout terms for 
leptons and sleptons in the last term in Eq. (J22i) . Correspondingly, for the density 
of the right-handed sneutrinos we use 

4^ 

^1,05 5^1,55 ^15 decay asymmetries for the decay of neutrino into Higgs 

and lepton, neutrino into Higgsino and slepton, sneutrino into Higgsino and lepton, and 
sneutrino into Higgs and slepton, respectively, defined by 

F~ — r~- F-.-._F~~ 

^i,a - sp _Lr~_^' ~ V fr ~ ~ + r ~ ~ 1 ■ ' ' 

In the MSSM, the four decay asymmetries are equal, e^ff^ = e^i^^ = ^sm _ ^mssm^ 
and given by [211 [25] 



.MSSM 



8vr (AU,)n \Mf 

with 



9 772 ' l^^ 



-In 



X \ X 



1 + x 



X 



(29) 



The matrix A depends on which MSSM interactions are in thermal equilibrium at the 
temperatures where leptogenesis takes place. Below (1 + tan^/5) x 10^ GeV, where the 
Boltzmann equations are solved for the individual asymmetries Fa^, and Ya^, A is 
given by 

-93/110 6/55 6/55 \ 
^MSSM ^ I 3/4Q _i9/3o 1/30 . (30) 

3/40 1/30 -19/30/ 

Between (1 + tan^/?) x 10^ and (1 + tan^/3) x 10^^ GeV, where the relevant flavour- 
dependent asymmetries are Ya2 = ^Ae+A^ and Fa^, we find 

.MSSM _ M41/761 152/761 \ 

~ V 46/761 -494/761 J ' ^"^^ 

and above (1 + tan^/5) x 10^^ GeV, we recover the flavour-independent treatment with 
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2.4 Solving the Boltzmann Equations in the SM and MSSM 



Solving the Boltzmann equations, for z from to oo, in the SM or in the MSSM yields 
the final B/?) — Loi asymmetries or 1^^^^ in the individual flavours. It is convenient 
to parameterize the produced asymmetries in terms of an efficiency factor Tja which, in 
the approximation that the small off-diagonal elements of A are neglected, is a function 
of AaaKa (no sum over a) and K, i.e. rja = ri{AaaKa, K), as 

Yl"^ = vT'el^^Y^-^izKKl), (32) 

^MSSM ^ ^MSSM^MSSM ^Y^^^^^<^^^^Yp^{z<^l)\, (33) 

generalizing the notation of [6] to the flavour-dependent treatment. Beyond the approx- 
imations of Sees. 12.21 and [2731 the efficiency factors also depend on Mi, and they depend 
on tan/3 in the MSSM. Y^{z -C 1) and Y^^z -C 1) are the number densities of the 
neutrino and sneutrino at T ^ Mi, if they were in thermal equilibrium, normalized 
with respect to the entropy density. In the Boltzmann approximation, they are given 
by 

Y^'liz « 1) - Yp^iz « 1) ^ -i^. (34) 



Eqs. ( 132|) and ( l33l) define the flavour- dependent efficiency factor in the SM and in the 
MSSM, and it can be computed by means of the Boltzmann equations fLOl [12] and its 
MSSM generalizations in Eqs. (1201) - (!22l) . where lepton flavour is taken into account 
correctly. The equilibrium number densities in Eqs. (1321) and (jHHj) serve as a normal- 
ization. A thermal population A^^i (and A'^i) decaying completely out of equilibrium 
(without washout effects) would lead to r^Q, = 1. Of course, K/Ka > 1 always holds 
by definition but K can be significantly larger than Ka- r]{AaaKa,K) is shown as a 
function of AaaKa in the MSSM for fixed values of K/\AaaKa\ = 2, 5 and 100 in Fig. [H 
In the SM, t] has the same qualitative (and a similar quantitative) behavior. 

What is more relevant than the differences in the flavour-dependent efficiency fac- 
tors (c.f. Fig. d]) is that the total baryon asymmetry is the sum of each individual lep- 
ton asymmetries, which is weighted by the corresponding efficiency factor [TUt [TT| [T^ . 
Therefore, upon summing over the lepton asymmetries, the total baryon number is 
generically not proportional to the sum over the CP asymmetries, ei = ei^a as 
in the flavour-independent approximation where the lepton flavour is neglected in the 
Boltzmann equations. In Eq. ([T^ . this corresponds to replacing Y^^ by the "one- 
single" flavour Ya = ^Ae+A^+A^ (the total lepton asymmetry), the flavour dependent 
decay asymmetries by ei, the washout parameters by K = Yla-^a^ matrix 
A by —1. The produced asymmetry in the flavour- independent approximation is then 
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Figure 1: Flavour-dependent efficiency factor ii{AaaKa, K) in the MSSM as a function of AaaKa, for 
fixed values of K/\AaaKa\ =2,5 and 100, obtained from solving the flavour-dependent Boltzmann 
equations [12 generalized to the MSSM (with tan/3 = 50, as an example), as displayed in Eqs. (PD|) 
- (HH). For larger K/\AaaKa\ the plot looks virtually like for K/\AaaKa\ = 100. More relevant 
than the differences in the flavour-dependent efficiency factors for different K/\AaaKa\ is that the 
total baryon asymmetry is the sum of each individual lepton asymmetries, which is weighted by the 
corresponding efficiency factor. As explained in the text, this can change the amount of the produced 
baryon asymmetry dramatically, compared to the flavour- independent approximation [101 1111 112j . 

given hj Y = eir]^'^'^{K)Y^{z <^ 1), where the flavour-independent efficiency factor 
r]^^'^{K) is related to ri{AaaKa, K) in Eq. ( l32l) by r7™'^(x) = r]{—x,x). In other words, 
in the flavour-independent approximation the total baryon asymmetry is a function of 
(Sa^i.") ^ ''^'"'^ (^/3-^/3)- correct flavour treatment the baryon asymmetry is a 

function of J2a^^,aV i^aaKa, K). If A''! decays about equally to all flavours and pro- 
duces about the equal asymmetry in all flavours, then (neglecting the effects of A for the 
moment) one expects that the flavour-independent approximation underestimates the 
asymmetry by roughly a factor of three for the case of strong washout and overestimates 
by a factor of three for the case of weak wash out. The approximate factor of three can 
be understood with the analytic approximations presented in \n\, from which we can 
see that in the case of weak (strong) washout for all flavours, the efficiencies are roughly 
(inverse) proportional to Ka- However, as we will see, there are situations in the SD 
models where the difference is much more dramatic. 

Let us note at this point how to generalize the above discussion for the range 
10^ GeV < Ml < 10^2 in the SM and (1 + tan^/?) x 10^ GeV < Mi < (1 + 
tan^ /3) X 10^2 GeV in the MSSM, where only the T is in equilibrium and is treated sep- 
arately in the Boltzmann equations, while the e and /i flavours are indistinguishable. In 
this case, following [I0l[12], one can combine the asymmetries for the e and /i flavours to 
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a combined density Y^^ = ^A^+Ap in the Boltzmann equations, where Ka is substituted 
by K2 = Ke + K^. The corresponding decay asymmetry is £1^2 = ^i,e + Above 
10^2 GeV in the SM and above (1 + tan2 13) x 10^^ GeV in the MSSM, we can combine all 
asymmetries for the e, /x and r flavours to a combined density Fa = FA<.+Ap+Ar in the 
Boltzmann equations and substitute Ka hj K = Ke + + K^-. The decay asymmetry 
then reduces to the flavour-independent one, Si = Si^e + + £i,t- 

The produced lepton asymmetries are partly converted into a final baryon asymmetry 
by sphalerons. For all temperature ranges, the produced baryon asymmetry (normalized 
to the entropy density) can be computed from the densities Y^ and Y^ssm gjyj- 
and MSSM, respectively, as 



yr = ^E^A^ (35) 
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3 Sequential Dominance 



To understand how sequential dominance works, we begin by writing the right-handed 
neutrino Majorana mass matrix Mrr in a diagonal basis as 

'Ma 

Mrr = I Mb I . (37) 



M, 



Cj 



Note that, as stated earlier, we work in a basis where the charged lepton Yukawa matrix 
is diagonal. In this basis we write the neutrino (Dirac) Yukawa matrix in terms of 
(1,3) column vectors Aj, _Bj, Cj as 

K={A B C) , (38) 

using left-right convention as in Eqs. ([T]) and ([2]). The Dirac neutrino mass matrix is 
then given by m^^^ = X^Vn- The term for the light neutrino masses in the effective 
Lagrangian (after electroweak symmetry breking), resulting from integrating out the 
massive right handed neutrinos, is 



_ {urA,){Aju,) {ufB,){Bju,) {ufC,){Cju,) 

Ma Mb ^ Mc ^ ^ 
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where z/j [i = 1,2,3) are the left-handed neutrino fields. Sequential dominance then 
corresponds to the third term being negligible, the second term subdominant and the 
first term dominant: 

Ma Mb Mc ' ^ ' 

In addition, we shall shortly see that small ^13 and almost maximal ^23 require that 

l^il < 1^2! ~ |^2|. (41) 

Without loss of generality, then, we shall label the dominant right-handed neutrino and 
Yukawa couplings as A, the subdominant ones as B, and the almost decoupled (sub- 
subdominant) ones as C. Note that the mass ordering of right-handed neutrinos is not 
yet specified. Again without loss of generality we shall order the right-handed neutrino 
masses as Mi < M2 < M3, and subsequently identify Ma, Mb, Mc with Mi, M2, M3 in 
all possible ways. The following results for the masses and mixing angles are independent 
of the mass ordering of the right-handed neutrinos. 

Writing = \Aa\e^'^^^ , = |-Ba|e"'^^i, Ca = \Ca\e^^^i and working in the mass 
basis of the charged leptons, under the SD condition Eq. ( HOl) . we obtain for the lepton 
mixing angles ^14j : 

tan ^23 ~ Yf \ , (42a) 

\B I 

tan ^12 ^ ^ , (42b) 

0231^2 I cos 02 - 5231^3 1 COS 03 

^ ,^ii^^..-^..^ \BiUlB, + AlB,) MA e-(^+^^.-^^.)|A| 
and for the masses 

^ (IA2P + IA3IV ..^ , 

ms ^ — , (43a) 

\B Pf^ 

mi ^ 0(|C|V/Mc) . (43c) 
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We would like to note that tri-bimaximal mixing [13] corresponds to the choice [TB] : 



l^ll 


= 0, 


\M 


= 1^43 


\Bi\ 


= \B2 




= 0. 



\B^ 



(44) 
(45) 
(46) 
(47) 



This is called constrained sequential dominance (CSD) [T6] . 

Let us now turn to the issue of leptonic CP phases in SD. In particular we will 
consider the MNS phase S, relevant for neutrino oscillations, and then study in Sec. 14.21 
how it is linked to CP violation required for leptogenesis. As in [19] the MNS phase S 
is fixed by the requirement that we have already imposed in Eq. fl42bp that tan(6'i2) is 
real and positive, 



C23|52|sin02 ~ S23I53I sin(^3 , 

C23I-B2I COS02 - S23I-B3I COS03 > 



(48) 
(49) 



where 



^B3 



+ 5. 



(50) 



The phase is fixed by the requirement (not yet imposed in Eq. fl42cl) ) that the angle ^13 
is real and positive. In general this condition is rather complicated since the expression 
for ^13 is a sum of two terms. However if, for example, Ai = then is fixed by: 

(51) 



(52) 



- 0Bi - C 

where 

C = aTg{A;B2 + A;B,). 
Eq. fl52|) may be expressed as 



tan(^ 



1^21 


\S23S2 + 


\B3\C23S3 


\B2 


S23C2 + 


-B3 C23C3 



(53) 



Inserting of Eq. flSTl) into Eqs. (jl8|) . (I50l) . we obtain a relation which can be expressed 

as 



tan(C + 6) 





^2 C23'52 — 


-33-52353 




\B2 


C23C2 + 


IBs 


\s23C3 



(54) 
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In Eqs. dSS]), dMD we have written Si = siiaQ, Ci = cosQ, where we have defined 

C2 = 0B2 - 0A2 , Cs = 0iJ3 - , (55) 

which are invariant under a charged lepton phase transformation. The reason why the 
seesaw parameters only involve two invariant phases rather than the usual six, is due 
to the SD assumption in Eq. ( l40l) that has the effect of effectively decoupling the right- 
handed neutrino of mass Mc from the seesaw mechanism, which removes three phases, 
together with the further assumption (in this case) of Ai = 0, which removes another 
phase. 



4 Flavour Matters in Leptogenesis with Sequential 
Dominance 

Let us now consider leptogenesis in neutrino mass models with sequential dominance 
(SD), taking into account lepton flavour in the Boltzmann equations. In SD, there are 
three classes of models with different characteristic predictions for leptogenesis. They 
differ by the role of the lightest right-handed neutrino in SD, which can either be the 
dominant one Mi = M^, the subdominant one Mi = Mb or the subsub dominant one 
Ml = Mc (which, in the SD limit, only contributes to mi but has a negligible effect on 
the neutrino mixing angles and CP phases). The possible form of the neutrino Yukawa 
matrix A^, is then given by 



A. 


= {A,B,C)otX,= 


{A,C,B), 


for Ml 


= Ma. 


A. 


= {B,A,C)otX,= 


iB,C,A), 


for Ml 


= Mb 


A. 


= {C,A,B)otX,= 


{C,B,A), 


for Ml 


= Mc. 



using the notation of Sec. [31 where we have ordered the columns according to Mrr = 
diag(Mi, M2, M3) where Mi < M2 < M3. The flavour specific decay asymmetries ei^a 
and washout parameters mi^„, calculated from Eqs. ( fT6l) and ( fTSll . are given in Tab. [1] 
for the three classes of SD. 



4.1 Decay Asymmetries and Washout Parameters 

We can now derive theoretical expectations and constraints on the flavour speciflc decay 
asymmetries and washout parameters from requiring consistency with the present neu- 
trino data in the classes of neutrino mass models discussed above. The experimental re- 
sults from neutrino oscillations indicate nearly maximal mixing 623 ~ 45°, large, but non- 
maximal mixing 612 ~ 33° and small mixing 613 < 13°, as well as small neutrino mass 
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Type of SD 






Ml = Ma 


3Mi 


Im 


AlB^iA^B)] ^ Im 


A^C^iA^C)] ] 


1 




Aa 


u 


167r 


Mb (At A) Mc (At A) 


Ml 


Ml = Mb 


3Mi 


'im 


'B*^A^iB^A)] ^ Im 


[B^C^iB^C)] 1 






Ba 




167r 


Ma {B^B) Mc {B^B) 


Ml 


Ml = Mc 


3Mi 1 




1 






2^2 


167r 


Ma{CW) Mb{C^C) 


Ml 



Table 1: Flavour specific decay asymmetries ei,a for the decay of the lightest right-handed neutrino 
with mass Mi (in the limit Mi <C M2,M3) and washout parameters toi.q in classes of models with 
sequential right-handed neutrino dominance (SD) in the SM. The MSSM decay asymmetries are a 
factor of 2 larger, in the considered limit. 

squared differences Am|i = m|— ^ 2.2x 10~^ eV^ and Am^i = ml — ml ^ 8.1x 10"^ 
eV^ [211 . In neutrino mass models satisfying SD, the neutrino Yukawa couplings and 
the masses of the right-handed neutrinos are related to the neutrino masses and mixings 
by Eqs. (iS]) and (iS]). 

We will focus first on the case where all three flavours are treated separately in the 
Boltzmann equations. As discussed in Sec. Ej this is the case for temperatures below 
10^ GeV in the SM and below (1 + tan^^g) x 10^ GeV in the MSSM. Typical examples 
where this case is relevant are thus unified models of flavour with large tan (3. We will 
then comment on the two remaining cases, where (i) the e and n flavour are combined 
in the Boltzmann equations and only the r is treated independently and (ii) where all 
flavours are combined and the treatment becomes flavour- independent, in Sec. 14.1.31 

4.1.1 Properties of the Decay Asymmetries ^i „ 

The decay asymmetries ei^a of Tab. [T] contain contributions from the two heavier right- 
handed neutrinos with masses M2 and M3. Using the SD conditions of Eq. ( 140|) . we 
see that the largest contributions to the decay asymmetries stem from the terms in 
Tab. [1] containing Ma, followed by the contributions containing Mb, and finally terms 
containing Mc give the smallest contributions which vanish in the limit mi 0. In 
the following, we will neglect the contributions to the decay asymmetry containing 
Mc, keeping in mind that they may be important when the leading contributions are 
suppressed, for instance due to a configuration of complex phases or due to a particular 
structure of the neutrino Yukawa couplings Ai and Bi. In the following, we will compare 
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the flavour specific decay asymmetries with the maximally possible decay asymmetry 

= , (59) 

which provides an upper bound on the decay asymmetry si in the SM in the flavour- 
independent approximation. In the MSSM, the maximal decay asymmetry is a factor of 
2 larger. The flavour specific decay asymmetries can be constrained in the three cases 
as follows: 



In the case Mi = Ma, assuming that the sub-dominant right-handed neutrino 
provides the leading contribution to the decay asymmetry and that the term con- 
taining Mc can be neglected, the flavour-specific asymmetries are given by 



^1. 



3Mi Im [AIB^{A^B)\ 
' 167r Mb {A^ A) 



(60) 



Using I All ^ 1 ~ I A3 1, as suggested by and Eqs. (I42cp and (142 al) and the 
smallness of ^13 together with the approximate maximality of ^23 , we can constrain 
Ei^p with (3 G r } and follows: 



< 



< 



3Mi Im [A*^Bp{AB)\ 



167r Mb (At A) 
3A'/il[|5J(|52| + 



167r 2 



m2 



o( 



Mb 
\Ai\ 
I A, I 



< 



15, 



3Mi [ 


^3| 


\B,\ 






1 + 


Bsl)] 


167r 




Mb2|A3| 


2 



\Bi\(\B2\ 



M 



B 



O 



( 



1^3 



O 



(61) 
(62) 



The estimate in the last step of Eqs. (16 ip and flU^ is based on the SD conditions 
in Eq. ( HOl) and on Eq. ( ]42bp . The latter states that in order to generate large 
mixing 6'i2 ~ 33°, |i?2| and/or |i?3| have to be of the same order as |-Bi|, which 
in turn is related to m2 by Eq. (143b[) . We conclude that |£:i_^| and |£:i,,-| are 
somewhat suppressed compared to the maximally possible value 5™**^ as was found 
for El = X]o^i,a in the flavour-independent treatment in [TH]. The value of |ei^e| 
given in Eq. flU^ is suppressed compared to |£i^^| and |ei^^| due to |Ai| -C IA2I, IA3I 
and it is related to the mixing angle ^13, as can be seen from Eq. (142 cl) . 
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Type of SD 


relation for ei g 


relation for ei^ 


relation for ei^r 


Ml = Ma 




< Q( '"a ymax 
~ ma ' 


< ^/ m2 ymax 
~ ma ' 


Ml = Mb 




^ _max 


_max 


Ml = Mc 




^ _max 


^ _max 



Table 2: Constraints and theoretical expectations on the flavour specific decay asymmetries from 
consistency with low energy neutrino observables in classes of models with SD in the SM, in the limit 
Ml <C M2, M3. e"^'^^, the upper bound on the decay asymmetry, is given by -j§:^^itT'3- In the MSSM, 
the decay asymmetries are a factor of 2 larger. 

• In the case Mi = Mb, using again |y4i| ^ IA2I ~ l^al, we find that we can 
constrain the decay asymmetries with /3 e {/i, t}, and ei^^ as 



3Mi Im [B;Ap{BU)] 



IGtt 



e - 



1^1. 



< 



Ma (BW) 
'\Bf,\{\B2\ + \B, 



< 



3Mi[|5^ 



\A,\\\B2\ + \Bs\ 



e - 



2 u. 



(BtE) 

|5i|(|52| + |53 



IGtt 



M^ (5tE) 



(fits) 



For l^i,/?!, the maximal value e™^^ can be nearly saturated, whereas |ei_, 
pressed due to \Ai \ <^ \A2\,\A^\ as in the case Mi = Ma- 



(63) 
(64) 

is sup- 



In the case Mi = Mc, the decay asymmetries ei^a obtain contributions from the 
right-handed neutrinos with masses Ma and Mb- Using analogous estimates than 
for the case Mi = Mb-, we can show that the contribution from the right-handed 
neutrino with mass Ma allows decay asymmetries |ei,^| and l^i.rl close to the up- 
per bound e™^^ similar to the result for Si in the flavour-independent treatment 
in [IH], but due to \Ai\ <^ IA2I ~ IA3I the decay asymmetry ei^e is suppressed. 
From Mb-, using estimates as for the case Mi = Ma-, all decay asymmetries can 
get contributions \si^a\ ^ 0{m2/m-!) e^^^. We will return to the case Mi = Mc in 
Sec. 14.3.31 where we discuss the specific example of tri-bimaximal mixing via con- 
strained sequential dominance (CSD). There we will also present explicit formulae 
for the decay asymmetries which illustrate the above statements. 



The above derived constraints on the decay asymmetries Si^a in the classes of models 
under consideration are summarized in Tab. [21 



18 



4.1.2 Properties of the Washout Parameters mi ^ 

• In the case Mi = M^, comparing the formula for in Eq. fl43ap with the formulae 
for the washout parameters mi ^ in Tab. [T] yields 

mi^^i = mi,T ~ . (65) 

Since for hierarchical light neutrino masses ms ~ ^/Am^ ^ 0.05 eV, we find 
that for these fiavours we are constrained to be in the strong washout regime, 
where mi ^ ^ m* (c.f Fig. [T]). On the contrary, m-i^e is significantly smaller 
since |y4i| <^ |742|, l^al and m-i^e is thus typically in the weak (optimal) region for 

mi e <^ m* (mi^e ~ m*). 

Note that in the flavour-independent approximation, the washout parameter would 
be given by mi = ~ ^s- The fact that we expect weak (or optimal) 

washout for rhg thus provides an important difference between flavour-dependent 
leptogenesis compared to the flavour-independent approximation. We will analyze 
this situation in more detail in Sec. 14.3.11 

• In the case Mi = Mb, from comparing the formula for m2 in Eq. fl43bl) with the 
formula for the washout parameter mi,e in Tab. [H, we infer that 

mi^e ~ Si2m2 . (66) 

Furthermore, from Eq. (142b[) and from the SD conditions Eq. ( l40l) it follows that 
B2 and/or S3 are of the order of Bi in order to generate large neutrino mixing 
^12, leading to the expectations 

mi^^ and/or mi^^ = 0(1712) . (67) 

• In the case Mi = Mc, Eq. (143 cl) implies that the washout parameters mi^a are 
typically < 0{mi) and therefore, with a hierarchical neutrino mass spectrum 
mi -C m2 < ma, they are in general in the weak (or optimal) washout region. 

A summary of the above derived constraints on the washout parameters mi^^ in the 
classes of models with sequential dominance is given in Tab. [31 
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Type of SD 


washout by rrii^e 


washout by mi^^ 


washout by mi 


Ml = Ma 


weak/optimal 


strong, fhi^^ ^ ^m^ 


strong, fhi^r ~ 


Ml = Mb 


optimal, fhi^e ~ Si2'^2 


< oim^r 


< o{m,y 


Ml = Mc 


weak/optimal 


weak/optimal 


weak/optimal 



Table 3: Flavour specific type of washout for the considered classes of models. (*At least one out of 
fni,ij.,fni,r has to be 0(7712) for Mi — Mb-) 

4.1.3 Remaining Cases: Only yT in Equilibrium and the Flavour-Independent 
Case 

As discussed in Sec. [21 all three flavours are only treated differently for temperatures 
smaller than 10^ GeV in the SM and larger than (1 + tan^ /?) x 10^ GeV in the MSSM. 

In the temperature ranges 10^ GeV < Mi < 10^^ GeV in the SM and (1 + tan^ P) x 
10^ GeV < Ml < (l + tan2/5) x lO^^GeV in the MSSM, only the yr is in equilibrium and 
thus only the r flavour is treated separately in the Boltzmann equations while the e and 
fi flavours are indistinguishable. In this temperature range, one can use the combined 
decay asymmetry and washout parameter 

£1,2 = £i,e + el^^, , (68) 
"^l,2 = "^i,e + rhi^^ , (69) 

instead of ei,e,£i,/x and mi,e,mi,/x and solve the Boltzmann equations for = ^a^+a^ 
and Fat- 

Above 10^2 in the SM and (1 + tan^/?) x lO^^ Q^y ^^e MSSM, all flavours 
are combined to the flavour-independent decay asymmetry and washout parameter 

£1 = £l,e + + £l,T , (70) 

fhi = fhi^e + fhi^^ + fhi^^ (71) 

and one can solve the Boltzmann equations for the total asymmetry Ya = ^Ae+A^+A^- 

The properties of the decay asymmetries and washout parameters in these two cases 
can be discussed analogously to Sees. 14.1.11 and 14.1.21 The results are presented in 
Tabs, a and El 
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Type of SD 


£1,2 




washout by mi_2 


washout by mi^^ 


Ml = Ma 




< Q( 'rri-2 ymax 


"ii,2 ~ l^^a 




Ml = Mb 






< 0{m,) 




Ml = Mc 


^ _max 


^ _max 


weak/optimal 


weak/optimal 



Table 4: Constraints and theoretical expectations on the flavour specific decay asymmetries and the 
flavour specific type of washout for the case where only the tau Yukawa coupling yr is in thermal 
equilibrium. For the e and /i flavours, the combined quantities £1^2 = £i,e + £i,fi and toi,2 = "^i.e + ^i./x 
are considered. 



Type of SD 


relation for 61 


washout by mi 


Ml = Ma 


< ^ ^ ^ ^max 


strong, mi = ms 


Ml = Mb 






Ml = Mc 




weak / optimal 



Table 5: Constraints and theoretical expectations on the flavour specific decay asymmetries and the 
flavour specific type of washout for the case where no charged lepton Yukawa coupling is in equilibrium. 
All flavours are combined effectively to a single flavour in the Boltzmann equations and the relevant 
quantities are the flavour-independent ones, Si — ei_e + + £i,t and rhi = fhi^e + "^-i,/^ + ™i.t- 

4.2 Leptonic CP Violation and the MNS - Leptogenesis Link 

In [19] it was shown how, in the case of SD where the lightest right-handed neutrino was 
the dominant one, there was a link between the MNS Dirac CP violating phase 6 which 
appears in neutrino oscillations and the phase which is relevant for leptogenesis. However 
the analysis in [19] was based on the flavour-independent formulation of leptogenesis, 
and it is therefore interesting to revisit this analysis in the light of flavour-dependent 
leptogenesis. 

Eqs. (1541) and fl33|) show explicitly that the MNS Dirac CP violating phase 5 is a 
function of only two seesaw phases (2, Csy ^is in [19]. Since these are the only remaining 
seesaw phases it must also be the case that leptogenesis must depend also on these 
same two phases. This opens up the possibility of a link between the leptogenesis phase 
and the oscillation MNS phase 6, as in [L9\- In the flavour-independent treatment of 
leptogenesis in [TH] the relevant leptogenesis phase was ( and the leptogenesis-MNS 
link was given by Eqs. (j53l) . (j54l) . However, in the present case, flavour-dependent 
leptogenesis will depend on both phases (2, Cs? we now discuss. 
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For the case Mi = M4, assuming that Ai = 0, we have 

^ 3M, lm[A*MA^B)] 
''"^ 167r Mb{AU) ■ ^ ' 

From this resuh it is clear that 
£^1,6 = 0, 

ei,^ oc +|A2||fi2| sin(C + C2), 

ei^r oc +|A3||S3|sin(C + C3). (73) 
From the previous section we see that 

Yb oc Fa, + Va. (74) 



with Ya given in Eq. flH^ . Since the efficiency factor rj^ depends on rhi^^ ^ 
\A2\'^v^/Ma and rj^- depends on rhi^r ~ {A^l'^v'^ / Ma, with |yl2| ~ {A^l we see that, 
r]f,^T]r, and 

oc 1^21 sin(C + C2) + 1^31 sin(C + Cs), (75) 

which shows exphcitly how leptogenesis depends on the two phases (2, Cs this 
case. 

For the case Mi = Mb, assuming that Ai = 0, we have 

^ 3M, lm[B:A^{B^A)] 
''""^ 167r Mb{BW) ■ ^ ' 

From this result it is clear that 

ei,e = 0, 

|A2||52|sin(C + C2), 

(77) 

Since now rj^ depends on the washout parameter mi,^ ^ \B2\^v^/Mb and r]^ de- 
pends on rhi j. ^ |i?3pf ^/M^, the two efficiencies are in general unequal in this 
case, and with IA2I ~ {A^l we see that, 

Yb oc -v^\B2\ sin(C + C2) - vAB^l sin(C + Cs), (78) 

which shows explicitly how leptogenesis depends on the two phases (2, Cs i^i this 
case. 



oc 


-1^2! 


1^21 


sin(C + C2 


£i,T oc 


-I^sl 




sin(C + Cs 
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• Clearly, for the case Mi = Mc, there is no leptogenesis-MNS link since leptogenesis 
will depend on phases associated with the right-handed neutrino of mass Mq, 
which does not contribute significantly to the seesaw mechanism, and therefore 
these phases will not significantly contribute to the MNS phases. 

The leptogenesis-MNS link in the flavour-dependent case is therefore rather similar 
to the link in the flavour-independent case, namely the connection between leptogenesis 
and 6 is made via the remaining two seesaw phases (2, Cs which are responsible for both 
leptogenesis and 6. These two phases fix 6 as in Eqs. (Ell), (!55l) . which are valid for all 
types of SD, including the cases Mi = Ma and Mi = Mb- However the two phases (2, Cs 
contribute to leptogenesis differently for the cases Mi = Ma and Mi = Mb, as shown in 
Eqs. ( ffSl) and ( !78l) . where both these results differ from the flavour- independent result 
in which leptogenesis is a simple function of (. 

In order to obtain a more precise leptogenesis-MNS link, one can reduce the number 
of phases still further by assuming an additional zero Yukawa coupling, in addition to 
assuming that Ai = 0. For example if we additionally assume that = 0, which 
implies that ei^r = 0, then this removes the phase (3, and we find that ( = (2- From 
Eqs. dlHl) - dSl, we obtain for the MNS phase S [12]: 

6 = -2C . (79) 

Under the same assumption that -B3 = (as well as Ai = 0), for the case Mi = Ma we 
see that 



oc +sin(2C), (80) 

while for the case Mi = Mb we see that 

rBCX-sin(2C), (81) 

where in both cases leptogenesis depends on a single phase (, which is directly related 
to the MNS phase 6 in Eq. ( I79l) . The sign of the CP violating phase 6 measurable 
in high precision neutrino oscillation experiments, or more precisely the sign of sin((5), 
therefore will depend on whether Mi = Ma or Mi = Mb- If Mi = Ma then sin((5) 
must be negative since baryon asymmetry of the universe in Eq. (180|1 must be positive. 
If Ml = Mb then sin((5) must be positive since baryon asymmetry of the universe in 
Eq. (IHTj) must be positive. Thus the measured sign of the oscillation phase 6 is capable 
of distinguishing between the two types of SD, namely Mi = Ma or Mi = Mb, under 
the assumption of the two texture zeroes Ai = and B3 = 00 



different possible choice of two texture zeros, which can be discussed analogously, is Ai ^ B2 = 0. 
In this case, the relation between 6 and C would be given by (5 = —2( + tt, leading to the prediction of 
the opposite sign of sin(5) from the observed baryon asymmetry. 
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4.3 Examples 



4.3.1 The Role of Ai in the Case Mi = Ma 

As we have remarked in the discussion of the flavour speciflc washout factors mi^aj 
one particularly interesting difference between flavour- dependent treatment and flavour- 
independent approximation is the possibility to have an optimal washout parameter mi e 
and strong washout in the flavours and r since in the flavour-independent treatment 
this leads to a strong total washout. We will therefore discuss the role of Ai in the case 
Ml = M^, sometimes also called light sequential dominance (LSD), in more detail. 

Let us start with the special case is Ai = 0, where due to |y42| ~ \A^\ required by 
6'23 ~ 45° (c.f. Eq. (142 al) ) the washout parameters satisfy 



iTLi^^ ~ Tni^r ~ '2^^ ' ■^^^'^^^S to rjfj, ^ rjr ^ 1] 



' A -— — 

. 2 m* ' m* 



12) 



We can now estimate for the baryon-to-entropy density in the MSSM, using Eqs. ([33 
and (EHD, 



10 / Img ms 

Yb ~ r] At-t , — 

31 V 2m* m* 



YZ{z » 1) 



r£(. » 1: 



(83) 



Apart from minor modiflcations in the efliciency factor, the flavour- dependent treatment 
is thus similar to the flavour- independent approximation in this case [T8] . 

The importance of a flavour-dependent treatment becomes clear if we allow for a non- 
zero Al <ti A2 ^ A3. As we have discussed above, rhi^e can now be optimal whereas mi,^ 
and fhi^r correspond to a strong washout. In the flavour-independent approximation, 
the total washout parameter rhi = Yla'^^,o^ ~ '"^s iniplies strong washout. In Fig. [21 
the flavour-dependent treatment is compared to the flavour-independent approximation 
in an example with non-zero Ai. Dashed lines correspond to the flavour- independent 
approximation and solid lines stand for the same examples with flavour-effects included. 
For non-zero A1/A2 above roughly 0.05, the produced baryon asymmetry is signiflcantly 
enhanced in the flavour-dependent case (more than two orders of magnitude in this 
example for A1/A2 = 0.2) because this induces an efliciency in the e-flavour which 
is close to optimal. In the flavour-independent approximation, the common efliciency 
factor 77™'^(mi/m*) is much lower due to the strong washout. 
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Figure 2: Example in the MSSM (with large tan/3 = 50), which illustrates the difference between 
flavour-dependent treatment and flavour-independent approximation for the case Mi = Ma if Ai is 
non-zero. We have chosen (j)Ai = —^,(j)B3 — Tr,4>A2 = = — 4'Ci = 0, -63 = 0,Ci ~ 0. The 
other parameters are determined from the experimental data assuming SD. For the parameters under 
consideration, ^13 remains below the present experimental upper bound. 



4.3.2 Leptogenesis and Tri-Bimaximal Mixing 

As discussed in Sec. El tri-bimaximal mixing can be realized in SD in a natural way 
via vacuum alignment and is sometimes referred to as constrained sequential dominance 
(CSD) [16J. As can be seen directly from Eq. (ST]), which states that A^B = B^A = 0, 
together with the formulae for the decay asymmetries ei^a (cf. Tab. [1]), in CSD the 
MSSM decay asymmetries reduce to 

For Ml = Ma and Mi = Mb, the decay asymmetries are thus suppressed by mi /ma. 

We would like to remark at this point that in realistic models, tri-bimaximal mixing is 
typically realized for the neutrino mass matrix in a basis where the charged lepton mass 
matrix is not entirely diagonal. However, including such charged lepton corrections do 
not change the above conclusions since A^B = B'^ A = remains unchanged by a change 
to the basis where the charged lepton mass matrix is diagonal, even though the lepton 
mixing is no longer exactly tri-bimaximal. In this respect, the conclusions are the same 
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in leptogenesis with correct flavour-dependent treatment and in the flavour-independent 
approximationj^ In general, however, we stress that charged lepton corrections do 
have an effect on leptogenesis in the flavour- dependent treatment, in contrast to the 
flavour- independent approximation. 

In contrast to the cases Mi = Ma and Mi = M^, tri-bimaximal mixing via CSD 
has excellent prospects for leptogenesis in the classes of models with Mi = Mc, as we 
now discuss. 



4.3.3 Tri-Bimaximal Mixing and CP phases in the Case Mi = Mc 

To be explicit, let us consider the speciflc choice A = (0, a, —a), B = {b,b,b) with real 
a and b and furthermore a texture zero in the (l,l)-element of A;^, i.e. C of the form 
C = (0, 6*2, C3). Then, ei e = and the remaining two decay asymmetries in the MSSM 
simplify to 

3Mi f Im[Qa^(Q-Q)] Im + Q)] ] 

Stt \ M^(CtC) MB(CtC) J 

3Mi f Im[Q(Q-Q)] ms Im[Q(Q + Q)] \ 

w 1 (cW) T + ioc) '^^""7 ' ^^^^ 

3Mi f lm[-C;a'iC;-C*,)] Im + Q)] ) 

Stt \ MA(CtC) MB(CtC) J 

3Mi [ Im[-C|(Q-Q)] m3 , Im[Q(Q + Q)] ^ ) 

w 1 pco ^ + {cW) '''"^'j ■ ^^^^ 

One can see that the two terms from the dominant (subdominant) RH neutrinos can 
easily contribute ~ ^max e™'^^s^2'^2/^3) to the decay asymmetries. This is nearly 
the optimal case for leptogenesis with hierarchical neutrinos. Furthermore, the washout 
parameters rhi^a are 0{mi), which allows optimal efficiencies rja if mi ^ m* ^ 10~^ eV. 

Let us now discuss under which conditions flavour effects are important, or play 
only a minor role. Similar to the discussion in Sec. 14.3.11 when the relevant washout 
parameters are equal, as this would in our example be the case for IC2I = IC3I, the 
situation is similar to the flavour-independent approximation. On the other hand, it 
is clear that if the washout parameters differ for the /i and r flavour, the produced 
asymmetry in the flavour- dependent treatment will differ in general from the flavour- 
independent approximation. 

^Note that in the R-matrix parameterization with R ~ 1, where aU mixings and phases in are 
introduced via the leptonic mixing matrix J/mnS; a-U flavour-specific decay asymmetries also vanish 
exactly as in the flavour- independent, since they are proportional to (Aj,A^)ij which is zero for 1 7^ J. 
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In addition, the values of the decay asymmetries depend strongly on the complex 
phases, in our simplified example on the two phases 0^2 (pa of C2 and C3, re- 
spectively. If IC2I 7^ IC3I and 0(^2 7^ 0C3, the differences between flavour-dependent 
treatment and flavour-independent approximation can be even more dramatic: To illus- 
trate this, let us focus on the dominant terms proportional to ms in Si^^ and Si^r given 
in Eqs. ( IHUj) and (!H7|) . which are proportional to 



£1,^ oc Im[C2*(C*-Q)], 
e,,^ oc \u^[-Cl{C; - CI)] . (89) 

On the other hand, in the flavour-independent approximation, the leading contribution 
to the decay asymmetry Si = ei^a = £i,fM + £i,t is proportional to 

E^ oc Im[(C*-Q)2]. (90) 

We can now imagine the situation that Im [(Cg — Cg)^] vanishes exactly, which would 
correspond to a choice of phases and (pca such that Arg(C2 — C3) G {0, |} (mod vr). 
Obviously, the phases (f)c2 &nd 0^3 can be chosen to satisfy this (for given IC2I 7^ IC3I) 
also with non-trivial 0^2 5 (pc'a ^ {0, f }, for instance with 



sin(0cj IC3 



(91) 



sin(0c3) IC2I 

which implies Im [(C| — C|)] = 0. Then, we have a situation where 

= -£i,T =■ £ (92) 

are non-zero, whereas the corresponding (leading) contribution to Si vanishes exactly. 
Since the washout parameters are different, the leading contribution to the decay asym- 
metries, using Eqs. (|33l) and (!35l) in the MSSM, can thus generates a baryon asymmetry 
of 

^ (r/, - Vr) i te(^ » 1) + rp^iz » 1)1 , (93) 
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proportional to the difference of the efficiencies r/^ = rj{A^^mi^^/m* , (^^ mi a)/?Ti*) and 
?7t- = ?7(y4T-T-mi^^/m*, (Xla^i.a)/"^*)- Let us consider the MSSM with tan/3 = 50 as an 
example, such that all flavours are to be treated independently for temperatures below 
about 2.5 X 10^^ GeV. The flavour-specific washout parameters fhi^^ and rhi^r are given 
by Clv^/Mi and Cffy/Mi, respectively (c.f. Tab. [I]). For instance, m^^ m* {Kr ~ 1) 
could lead to a nearly optimal efficiency factor 1]^-, while C2 <^ C3, and correspondingly 
f^* {K^ <^ 1), could imply a very low efficiency 1]^ <^ rj^ (c.f. Fig. [T]). According 
to Eq. fl93p . in the correct flavour-dependent treatment a baryon asymmetry can be 
generated, while in the flavour-independent approximation (in leading order in SD) it 
would vanish. This situation is illustrated in Fig. [3l 
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Figure 3: Evolution of the flavour specific B/3—La asymmetries Ya^ and Ia^ as a function oi z — Mi/T 
for the example of tri-bimaximal mixing via CSD, with Mi = Mc and the choice Ci = 0. Complex 
phases cjyc^ and (/iCg are chosen as in Eq. leading to ei,r — ~£i,n- While the flavour-independent 
approximation yields zero lepton asymmetry Find, the two flavour-specific asymmetries evolve differently 
for different flavour-specific washout parameters = mi^^/m* and K^- = rhx^rj'fn* ^ allowing a final 
baryon asymmetry to be produced (c.f. Eq. (|93p ). This situation can be realized, for example, with 
Ml = Mc = 1.85x10" GeV, |Ci| = 0, IC2I = 1.69x10-4, IC3I = 3.09 x lO^^, = -7r/2,0C3 = 0.028 
and Ma = 5 X lO^^ GeV, |^i| = ^,\Ai\ = |^3| = 0.063, (/ia^ = ^Aa^ = 0. The subdominant SD 
contribution to the decay asymmetries vanishes for = 0,(13 B3 = ^Ca- The decay asymmetries are 
given by ei_T ~ —10^^, ei ^ 10^^, and the washout parameters are w 1, Kf^ w 3 x 10^"^. 

5 Summary and Conclusions 

We have studied thermal leptogenesis in a class of neutrino mass models based on the see- 
saw mechanism and SD, taking into account lepton flavour in the Boltzmann equations 
in both the SM and MSSM. SD models are a very well motivated class of models which 
have been widely applied to unified flavour models, and thus are quite representative 
of the general class of three right-handed neutrino models where the importance of 
flavour- dependent leptogenesis effects has hitherto not been evaluated. 

In general, flavour-dependent effects of leptogenesis are relevant at low enough tem- 
peratures (set by the lightest right-handed neutrino mass) such that at least one charged 
lepton flavour is in thermal equilibrium. When this condition is met, flavour-dependent 
effects are important when the efficiency factors rja differ significantly for the distin- 
guishable flavours. The efficiency factors rja depend on rn-i „ which in SD are determined 
from the magnitudes of the Yukawa couplings of the lightest right-handed neutrino to 
the different lepton flavours, as indicated in Tab. [H 
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If the dominant right-handed neutrino is the hghtest one, then only Tjg may differ, 
since \A2\ ~ {A^l imphes that r)n ^ rjT, and flavour-dependent effects are only important 
when Ai ^ 0, corresponding to non-zero texture element in the (l,l)-element of the 
neutrino Yukawa matrix X,y. We have studied this case numerically and shown that 
dramatic effects can result. Moreover, we have also seen that, even with Ai = 0, 
flavour-dependent effects are relevant for the leptogenesis-MNS link since the decay 
asymmetries ei_^, ei^T may differ due to B2 ^ B^. An extreme example of this is when 
there are two zero textures in the dominant and leading subdominant parts of X^,, for 
instance = B3 = 0, in which case the cosmological CP violating phase for leptogenesis 
is directly related to the CP violating phase 5 of the MNS matrix observable in neutrino 
oscillations. Furthermore, the baryon asymmetry of the universe under this conditions 
determines the sign of the MNS phase 5 (i.e. the sign of sin(5)) for the case where the 
hghtest right-handed neutrino is the dominant one. We have shown that this phase is 
opposite to the case where the lightest right-handed neutrino is the leading subdominant 
one, enabling the two types of SD to be distinguished experimentally. 

For the case that the lightest right-handed neutrino is the almost decoupled one, there 
is no leptogenesis-MNS link. Nevertheless this case is amongst the most promising for 
thermal leptogenesis since the efficiencies may be optimal. Although this is also true 
in the flavour- independent treatment, flavour-dependent effects are expected to play an 
important role when the efficiences are unequal, corresponding to the magnitudes of the 
Yukawa couplings \Ca\ being unequal. Unlike the previous case, this can be achieved 
with a zero texture in the (l,l)-element of Aj, corresponding to Ci = 0, since here we can 
have IC2I ^ IC3I since these couplings are unconstrained by neutrino phenomenology, 
and so in general may be quite different. Thus very large flavour-dependent effects 
can be present in this case, an extreme example being when ei.^ and ei.^ arc equal 
and opposite which would lead to a zero result in the flavour-independent case, but a 
non-zero result in the flavour-dependent case due to the different efficiency factors. 

In conclusion, we have seen that in many cases ffavour- dependent effects may be 
important for leptogenesis in SD models involving three right-handed neutrinos due 
to the efficiencies of the distinguishable flavours being different. The effects range from 
factors of three, up to differences of a few orders of magnitude, or in extreme cases having 
a large non-zero ffavour-dcpcndent result for cases where the ffavour-independent result 
gives zero. In other cases the flavour-dependence of leptogenesis is not important for the 
asymmetry, due to the efficiencies of distinguishable flavours being equal, as in the case 
when the lightest right-handed neutrino dominates the see-saw mechanism, and there 
is a zero texture in the (1,1) -element of X^,. However even in this case the relationship 
between leptogenesis phases and the MNS phase is affected by flavour-dependent effects 
due to the lepton asymmetries being unequal. We conclude that flavour-dependent 
effects cannot be ignored when dealing with three right-handed neutrino models. 
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